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1. Introduction 

The fields of statistical physics, condensed matter physics and quantum information 
theory share a common interest in the study of interacting quantum many body systems. 
The concept of entanglement in quantum mechanics has significant importance in all 
these areas. Much of the current effort is devoted to the description and quantification 
of the entanglement contained in strongly correlated quantum states. Quantum 
entanglement is a fundamental measure of how much quantum effects we can observe and 
use to control one quantum system by another, and it is the primary resource in quantum 
computation and quantum information processing [9], [61]. Entanglement properties 
play an important role in condensed matter physics, such as phase transitions |68j . 
[69] and macroscopic properties of solids [30]. Extensive research has been undertaken 
to investigate quantum entanglement for spin chains, correlated electrons, interacting 
bosons as well as other models, see [3], [6], [8], [H], [H], [H], [31], [32], [39], [12], [13], 

isu, inzi, m 

[81], [HI] for reviews and references. Characteristic functions of quantum entanglement, 
such as von Neumann entropy and Renyi entropy, were obtained and discussed through 
studying reduced density matrices of subsystems [19], [26], [27], [16], [1^. An area law 
for the von Neumann entropy in harmonic lattice systems has been extensively studied 



Much insight in understanding entanglement of quantum systems has been obtained 
by studying exactly solvable models in statistical mechanics. In 1987, I. Affleck, T. 
Kennedy, E. H. Lieb and H. Tasaki proposed a spin interacting model known as the 
AKLT model [1], [2]. The model consists of spins on a lattice and the Hamiltonian 
describes interactions between nearest neighbors. The Hamiltonian density is a linear 
combination of projectors. The model is similar to the Heisenberg anti-ferromagnet 
with a gap. The authors (AKLT) of [1], [2] found the exact ground state, which has 
an exponentially-decaying correlation function and a finite energy gap. This model has 
been attracting enormous research interests since then [I^, [17], [21], [50], [45]. It can 
be defined and solved in higher dimensional and arbitrary lattices [2] , [20] , [52] , [71] and 
generalizable to the inhomogeneous (non-translational invariant) case (spins at different 
lattice sites may take different values) and an arbitrary graph [53]. Given certain 
conditions (as to be described later), the ground state has proven to be unique [5], [53] . 
It is known as the Valence-Bond-Solid (VBS) state. The Schwinger boson representation 
of the VBS state (see (fT7|l ) relates to the Laughlin ansatz of the fractional quantum Hall 
effect [5], [37], [39], [H], [59]. The Laughlin wave function of the fractional quantum Hall 
effect is the VBS state on the complete graph [12], [36]. The VBS state illustrates ground 
state properties of anti-ferromagnetic integer-spin chains with a Haldane gap [35] . The 
theory of VBS state was essentially developed by B. Nachtergaele and others [22], [23] , 
[21], [22], [63], [M]- The Entanglement of formation in VBS state was estimated in 
[62] . Brennen and Miyake showed that the VBS state can be used as a resource state in 
measurement-based quantum computing instead of the cluster state plij. It was proved 
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in that VBS state allows universal quantum computation and an implementation 
of the AKLT Hamiltonian in optical lattices |29] has also been proposed. 

We shall consider a part of the system, i.e. a block of spins. It is described by the 
density matrix of the block, which we call the density matrix later for short. The density 
matrix has been studied extensively in [19], [28], [19], [53], [77]. It contains information 
of all correlation functions [5], [17|, [19], [82]. Furthermore, the entanglement properties 
of the VBS ground state has been studied by means of the density matrix as in [T7] . 
|19] . [20] . |42j . [19] . [77] . The von Neumann entropy of the subsystem density matrix is 
a measure of entanglement of the VBS state. The Renyi entropy is another measure of 
the entanglement. The entanglement entropy was obtained in [19], [28], [19], |82j . 

The structure of the density matrix is important. For a 1-dimensional AKLT spin 
chain the density matrix has a lot of zero eigenvalues [82] , [83] . The eigenvectors with 
non-zero eigenvalues are the degenerate ground states of some Hamiltonian, which we 
shall call the block Hamiltonian (see (123!) ). In the limit of large block, the density matrix 
is proportional to a projector on the degenerate ground states of the block Hamiltonian. 
These states are the only eigenstates of the density matrix with non-zero eigenvalues 
which contribute to the entropy. Also, eigenstates of the density matrix are useful in 
quantum computing because of quantum measurements. It was conjectured in [82] that 
eigenvectors of the density matrix with non-zero eigenvalues always form degenerate 
ground states of some Hamiltonian (the block Hamiltonian), which is generalizable to 
an arbitrary graph. In this paper we shall give a general proof of this statement. 

The paper is divided into five parts: 

(i) We define the basic AKLT model on an arbitrary connected graph and construct the 
unique VBS ground state using symmetrization and anti-symmetrization of states. 
A graphical illustration is included. (Section [2]) 

(ii) We introduce the generalized (inhomogeneous) AKLT model and give the condition 
of the uniqueness of the ground state. The VBS ground state is constructed using 
the Schwinger boson representations. Within this formulation, the relation between 
the VBS state and the Laughlin states of the factional quantum Hall effect becomes 
obvious [5], [12], [36], [37], [39], [59]. (Section[3D 

(iii) In order to study entanglement of the VBS state, we cut the graph into 
two disconnected parts: the block and the environment. We define the block 
Hamiltonian, and show that its ground state is degenerate. (Section H]) 

(iv) The density matrix of the block is proved to have a lot of zero eigenvalues. The 
eigenvectors with non-zero eigenvalues form degenerate ground states of the block 
Hamiltonian. (Section [5]) 

(v) Examples of the density matrix are given explicitly as special cases of the general 
result. We also formulate some open problems. (Section [6]) 

In Sections [2] and [3] we follow the paper of A. N. Kirillov and V. E. Korepin written 
in 1990, see [53]. 
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2. The Basic AKLT Model 

We start with the definition of the basic AKLT model on a connected graph. A graph 
consists of two types of elements, namely vertices and edges. Every edge connects two 
vertices. As in Figure ([T]), a vertex is drawn as a (large) circle O and an edge is drawn 

as a solid line connecting two vertices. For every pair of vertices in the connected 

graph, there is a walk from one to the other. Vertices can also be called sites and edges 
sometimes called links or bonds. In the case of a disconnected graph, the Hamiltonian 
([1]) is a direct sum with respect to connected components and the ground state is a direct 
product. We shall start with a connected graph. We shall also assume that the graph 
consists of more than one vertices, otherwise there would be no interaction at all. Let's 
introduce notations. By Si we shall denote the spin operator located at vertex / with 
spin value Si. In the basic model we require that Si = \zi, where zi is the number of 
incident edges (connected to vertex /), also known as the valence or coordination number 
(the number of nearest neighbors of the vertex /). The relation between the spin value 
and coordination number must be true for any vertex /, including boundaries. This 
will guarantee uniqueness of the ground state. The Hamiltonian describes interactions 
between nearest neighbors: 

H = J2H{k,l)- (1) 

(kl) 

Here H{k, I) describes the interaction between spins at vertices k and I connected by 
an edge, and we sum over all edges {kl). The Hamiltonian density is H{k,l). To write 
down an explicit form of H{k,l), we define a projector 7ij{k,l): 

Ik n TT (s> + s,)' - id + 1) 

Operator 7ij{k,l) projects the edge spin Jm = Sk + Si on the subspace with fixed 
total spin value J and |S'a; — < J < Sk + Si. Note that we could expand 
{Sk + Si^ = 2Sk ■ Si + SkiSk + 1) + Si{Si + 1). So that projector 7rj{kJ) in is 
a polynomial in the scalar product {Sk ■ Si) of degree 2Smin, where Smin = Min{Sk, Si} 
is the minimum of the two spin values of the same edge. For example with Sk = Si = 1, 
we may have a quadratic polynomial: 

MkJ) = l{Sk-Siy + ^{Sk-Si) + ^. (3) 

In the basic model we define the Hamiltonian density H{k, I) as 

H{k, I) = A{k, I) 7vs,+s, {k, I), H{k, I) > (4) 

with A{k,l) an arbitrary positive real coefficient (it may depend on the edge (kl)). 
So that the Hamiltonian density for each edge is proportional to the projector on the 

I A walk is an alternating sequence of vertices and edges, beginning and ending with a vertex, in which 
each vertex is incident to the two edges that precede and follow it in the sequence, and the vertices 
that precede and follow an edge are the endvertices of that edge. 
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Figure 1. Example of a part of the graph including vertex k with ~ 3 and vertex I 
with zi = 4. Black dots • represent spin-^ states, which are enclosed by large circles 
O representing vertices and symmetrization of the product of spin-^'s at each vertex. 
Solid lines represent edges which anti-symmetrize the pair of connected spin-^'s. 



subspace with the highest possible edge spin value {Sk + Si). The physical meaning is 
that interacting spins do not form the highest possible edge spin (this will increase the 
energy) in the ground state. The Hamiltonian in ([1]) is a linear combination of projectors 
with positive coefficients, which shows that H is semi-positive definite. 
The Hamiltonian ([T]) with condition 

Si = (5) 

has a unique ground state [1], [2], [5], [53] known as the Valence-Bond-Solid (VBS) state. 
It can be constructed as follows. Each vertex / has zi spin-^'s. We associate each spin-^ 
with an incident edge. In such a way each edge has two spin-|'s at its ends. We anti- 
symmetrize the wave function of these two spin-|'s. So that anti-symmetrization is done 
along each edge. We also symmetrize the product of spin-^'s at each vertex (each large 
circle). Let's write down the VBS ground state algebraically. We label the particular 
dot from vertex I connected with some dot from vertex k by Ik (correspondingly, that 
dot from vertex k is labeled by ki). In this way we have specified a unique prescription 
of labels with dots. Then the anti-symmetrization results in the singlet state 

mki = ^{\ i)k,-\iW\ T)fcJ. (6) 
The direct product of all these |$) singlet states corresponds to all edges in our graph: 

(7) 

{kl) 

We still have to complete the symmetrization (circles) at each vertex. We denote the 
symmetrization operator of Zi dots in vertex / by P(05 then the symmetrization at each 
vertex is carried out by taking the product n«P(0 of all vertices. Finally, the unique 
VBS ground state can be written as 

|VBS) = ^P(O^I*)'c^ (8) 

/ {kl) 



VBS on Arbitrary Graph 



6 



Here the first product runs over all vertices and the second over all edges. If the 
coordination number zi is a constant over all vertices in the graph except for boundaries, 
then we would have the same spin value at each bulk vertex. In that case the basic model 
is also referred to as the homogeneous model. 



3. The Generalized AKLT Model 



In the generalized AKLT model, relation ([5]) is generalized. We associate a positive 
integer M^i {M^i = Mi^) to each edge {kl) of the graph. We shall call M^i multiplicity 
numbers. The Hamiltonian describes interactions between nearest neighbors (vertices 
connected by an edge): 

H = J2H{k,l)- (9) 

(kl) 

However, the Hamiltonian density is no longer proportional to a single projector in 
general. It is a linear combination of projectors 

S k+S i 

H{k,l)= Yl Aj{k,l)nj{k,l), H{k,l)>0. (10) 

J=Sk+Si-Mki+l 

Projector 7rj{k,l) is given by ([2]), and Aj{k,l)'s are arbitrary positive coefficients. So 
that H{k, I) projects the edge spin on the subspace with spin value J greater than 
Sk + Si — Mki. Physically formation of edge spin higher than Sk + Si — M^i would 
increase the energy. Cappelli, Trugenberger and Zemba showed that the Hamiltonian 
for the fractional quantum Hall effect also can be written as a linear combination of 
projectors [T2] . 

The condition of uniqueness of the ground state was introduced in [53]: 

2Si = J2Mki, V/. (11) 

k 

Here Si is the spin value at vertex / and we sum over all edges incident to vertex / 
(connected to vertex I). The Hamiltonian ([9]) has a unique ground state if (|TT1) is 
valid. The relation 25"; = zi for the basic model is a special case when M^i = 1. The 
condition ( ITTi) can be put into an invariant form. Let's define a column vector S, the l^^ 
component of which is associated with vertex / of the graph and equal to Si. The number 
of components is equal to the number of vertices N. Next, we define another column 
vector M with its dimension equal to the number of edges M in the graph. The k^^ 
and components of this vector are associated with edge (kl) and both equal to M^.;. 
The most important geometrical characteristic of the graph is the vertex-edge incidence 
matrix I (see |40]|). This is a rectangular matrix with rows and M columns. Each 
row is associated with the vertex and each column is associated with the edge. If the 
vertex belongs to the edge the corresponding matrix element is equal to one, otherwise 
zero. Then the condition ffTTl) of uniqueness can be re-written as 

2S = i-M. (12) 

For more details we refer to ^53j . 
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Under condition (ITT]) or (1121] . the unique ground state of Hamiltonian (Q is referred 
to as the generahzed VBS state. It is constructed by introducing the Schwinger boson 
representation [5], [28], [53], [19], [82], [H3]. We define a pair of independent canonical 
bosonic operators a/ and bi for each vertex /: 

[ ttfc, a] ] = [ bk, b\ ] = 6ki (13) 
with all other commutators vanishing: 

[ ttfc, ai] = [ bk, bi] = [ak, ] = [ Qfc, b] ] = 0, V A;, I. (14) 
Spin operators are represented as 

S+ = at6,, Sr = bjau Sl = ^iajai-b}k). (15) 

To reproduce the dimension of the spin-S"/ Hilbert space at vertex /, a constraint on the 
total boson occupation number is required: 

^{alai + b}k) = Si. (16) 
As a result, the VBS ground state in the Schwinger representation is given by 

\VBS) = l[{4bl-biaiy'''\ye.c), (17) 

(kl) 

This representation shows that for a full graph (each vertex is connected to every other 
vertex) the VBS state coincides with the Laughlin wave function [S], [12], [3S], [SZ], [3S]- 
In f fT7|) the product runs over all edges and the vacuum |vac) is annihilated by any of 
the annihilation operators: 

ai |vac) = bi |vac) = 0, V /. (18) 

Note that [ al, b] ] = 0, W k,l. If we replace a^'s and 6^'s by complex numbers 
(coordinates of electrons), then (ITTll turns into the Laughlin wave function of the 
fractional quantum Hall effect [5], [37]. To prove that f|T7j) is the ground state we 
need only to verify for any vertex / and edge {kl): (i) the total power of a] and b] 
is 25*;, so that we have spin-S*/ at vertex /; (ii) — |(I]//^/M//fc + J2k'^k Mk'i) < J^i = 
^k + < ^(J2i'^i Milk + J2k'^k ^k'l) by a binomial expansion, so that the maximum 
value of the edge spin Jki is |(E/vi M'fc + Y^k'^kMyi) = Sk + Si - Mki (from SU{2) 
invariance, see 0). Therefore, the state |VBS) defined in f fT7|) has spin-S'^ at vertex / 
and no projection onto the Jki > Sk + Si — Mki subspace for any edge. The introduction 
of Schwinger bosons can be used to construct a spin coherent state basis in which spins 
at each vertex behave as classical unit vectors. The spin coherent state basis is given in 



Appendix A 



4. The Entanglement between Block and Environment 

The VBS state (see and ( ITTjl ) has non-trivial entanglement properties. The density 
matrix of the VBS state is a projector 
|VBS)(VBS| 



(VBS|VBS) 



(19) 
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Let us cut the original graph into two subgraphs B and E, that is, we cut through 
some number of edges such that the resuhing graph B U E is disconnected (no edge 
between B and E). We may call one of them, say B, the block, and the other one E 
the environment. The distinction is arbitrary and the two subsystems are equivalent in 
measuring entanglement. 

Let's focus on the block (subsystem B). It is described by the density matrix of 
the block (obtained by tracing out degrees of freedom of the environment E from the 
density matrix p ( fT9i) ): 

Pfe = Tre[p]. (20) 
In fl20|) and below we use subscript b for block and e for environment. The density 
matrix ph contains all correlation functions in the VBS ground state as matrix entries 
0, [S], [l9], [82]. The entanglement can be measured by the von Neumann entropy 

S'^.AT = -Tr4 Pfelnpb ] = - ^ AlnA (21) 
or the Renyi entropy 




Here A's are (non-zero) eigenvalues of density matrix ph and a is an arbitrary parameter. 
It was shown by using the Schmidt decomposition [65j that non-zero eigenvalues of 
the density matrix of subsystem B (block) is equal to those of the density matrix 
of subsystem E (environment). So the two subsystems are equivalent in measuring 
entanglement in terms of entanglement entropies, i.e. Sy,]\f[B] = Sy,]y[E] and Sji[B] = 
Spi[E]. This fact has been used in obtaining entanglement entropies of 1-dimensional 
VBS states as in [19], [19]. We shall show that the spectrum of the density matrix pb 
contains a lot of zero eigenvalues. In order to understand this and give the subsystem 
(block) a more complete description, we first introduce the Hamiltonian of the subsystem 
(called the block Hamiltonian). 

The block Hamiltonian Hb is the sum of Hamiltonian densities H{k, I) with both 
k E B and / G B, i.e. nearest neighbor interactions (edge terms) within the block B: 

Hb= H{k,l), keB, leB. (23) 

Here H{k, I) is given in (jl]) for the basic model and (ITO!) for the generalized model, 
for k and / connected by an edge. In (|23l) no cut edges are present (boundary edges 
between subgraphs B and E removed). This Hamiltonian has degenerate ground states 
because uniqueness conditions (I5l) and ( ITTi) are not valid. Let us discuss the degeneracy 
of ground states of (l23l) . Let's denote by L the number of vertices on the boundary of 
the block B. The boundary consists of those vertices with one or more cut incident edge, 
see Figure ([2]). The degeneracy deg. of ground states of Hb is given by the Katsura's 
formula [48j 



deg. = Y[ 

ledB 



J2 Mki]+l 



{kl) e {cut edges}. (24) 
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Figure 2. Example of the cutting for the basic model. The curved double line 
represents the boundary between the two subgraphs. We have the block B on the left 

and the environment E on the right. Solid lines represent edges while dashed 

lines represent cut edges. Each dashed line connects two dots. All vertices in 

the figure belong to the boundary of B or E because of the presence of one or more 
cut incident edges (dashed lines). 



Here dB denotes vertices on the boundary of the block B and dE are vertices on the 
boundary of the environment E. In ( 124|) we have L terms in the product. Formula 
(!24|) is valid for both the basic and the generalized model. For the basic model all 
Mki = 1. A proof of formula ( !24l) is given in Appendix B The subspace spanned by 
the degenerate ground states is called the ground space, with the dimension given by 
deg. in flM|) . We emphasize at this point that the block B should contain more than 
one vertices, otherwise the block Hamiltonian vanishes Hb = and the whole Hilbert 
space become the ground space. We discuss the density matrix for a single vertex block 
at the end of next section. It was shown for 1-dimensional models in [82], [M] that the 
spectrum of density matrix ph is closely related to the block Hamiltonian. The density 
matrix is a projector onto the ground space multiplied by another matrix. We shall 
prove the statement for an arbitrary graph in the next section. 



5. The Density Matrix 

Let us denote by A*";, the number of vertices in the block B. Then the dimension dim. 
of the Hilbert space of the block B is equal to Yli{2Si + 1) with / G B, which is also the 
dimension of the density matrix p^. The value is 

dim. = l[[zi + l], (25) 

l£B 

for the basic model and 
dim. = 

leB 



Mm\+1 



vfee(-Buas) 



(26) 
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for the generalized model. In both expressions fl25p and (1261) we have A*";, factors in the 
product. The density matrix would have dim. number of eigenvalues. However, most 
of the eigenvalues are vanishing and pb is a projector onto a much smaller subspace 
multiplied by another matrix. To prove the statement, we define a support to be 
the subspace of the Hilbert space of the block B with non-zero eigenvalues, i.e. it 
is spanned by eigenstates of pb with non-zero eigenvalues. The dimension of the support 
is denoted by D. We have the following theorem on the structure of the density matrix 
Pb (Assuming that the block have more than one vertices, i.e. Nb > 2, so that Hb is not 
equal to zero identically): 

Theorem The support of pb (l20l) is a subspace of the ground space of the block 
Hamiltonian Hb fl23l) . 

To prove the theorem, we recall that H = Yj{ki)eB H{k, I) and each H{k, I) is a sum 
of projectors ( |TOl) . We have H{k,l) > 0. Then the construction of the VBS ground 
state ([8]) and (fT7|) guarantees that there is no projection onto the subspace with higher 
edge spins for any edge. Therefore, 

H{k,l)\VBS) =0, \/ (kl). (27) 

In particular, this is true for edges inside the block B, i.e. k E B and / G B. Now, from 
the definition of pb in ( 120|) . we have 

H{k,l)pb = H{kJ)TT,[p] 

_ H{kJ)TTe[ |VBS)(VBS| ] 
~ (VBS|VBS) 

In the last step of (125]) we have used ([27]) and the fact that edge (kl) lies completely 
inside block B so that H{k, I) commutes with tracing in the environment E. Equation 
f l28|l is true for any edge in 5, so that 

HbPb= E H{k,l)pb = 0, keB, lEB. (29) 

(kl)£B 

If we diagonalize the density matrix pb 

p, = EA|A)(A|, (30) 

where |A) is the eigenstate corresponding to eigenvalue A. Then (!29|) can be re- written 

as 

i/,5:A|A)(A| = 5:Ai/,|A)(A|=0, (31) 

A^O At^O 

Note that {|A)} is a linearly independent set. Therefore the solution of (13T|) means that 

Hb\X) = 0, X^O. (32) 

Expression (152]) states that any eigenstate of pb (with non-zero eigenvalue) is a ground 
state of Hb- As a result, we have proved the Theorem that the support of pb is a 
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subspace of the ground space of Hh, so that D < deg. The density matrix takes the form 
of a projector multiphed by another matrix and the projector projects on the ground 
space. Also, it is clear from expression ( 12^ and f l25i) . ( l26l) that deg. < dim. {dB C B 
so that L < Nh). Usually, deg. is much smaller than dim. because the former involves 
only contributions from boundary vertices of the block while the latter also involves 
contributions from all bulk vertices. Then as a corollary of the Theorem, we have 
D < deg. < dim. 

If the block B consists of only one vertex with a spin-S", then we conjecture that it 
is in the maximally entangled state. The support has dimension D = 2S + 1. 



6. Examples of the Density Matrix and Open Problems 



The density matrix of the block has been studied in |19], [49] and diagonalized directly 
in [82], [83] for 1-dimensional models, which illustrates the Theorem explicitly. It was 
shown for different 1-dimensional AKLT models that the inequality D < deg. is always 
saturated, i.e. D = deg., so that the support is exactly equal to the ground space. The 
density matrix is proportional to the projector on the degenerate ground space of the 
block Hamiltonian. Therefore the projector P^ on the support of pb is equal to the 
projector Pdeg. on the ground space of Hi,: 

Pd = Pdeg. (33) 

If we denote the identity of the Hilbert space of the block by Idim., then we also have 

Pbihim. - Pdeg.) = 0. (34) 

Using these relations fl33l) and flMl) . the density matrix can be put in the following 
matrix form 

Pf, = A • Pdeg. = Pdeg. ' A (35) 

where A is a diagonal matrix with non-zero eigenvalues of pt, as entries. It was also 
proved in [82], [83] that in the large block limit Ni, — ^ oo, all eigenvalues become the 
same so that 

hm A=^ln, (36) 

where Id is the identity of the support. As a consequence, the density matrix approaches 
the following limit 

lim pb = Poo = -y—Pdeg. = ^Pd, (37) 
Nt^oo' deg. D 

where poo behaves as the identity in the ground space or support. 

In below we give explicitly the forms of the density matrix obtained in [82], [83] 
and formulate some open problems. 
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6.1. 1- dimensional Basic Model 



For the basic model in 1-dimension, we have spin-l's in the bulk and spin-|'s at both 
ends of the chain. The block B consists of contiguous bulk vertices and the block 
Hamiltonian is 

(Af(,-1) terms 



(3^ 



E 

(I, l+l) e block 

As shown in p2j, the ground space of Hh is 4-dimensional. It can be spanned by 
{\G]a),a = 0,1,2,3}, and these are also eigenstates of pb with non-zero eigenvalues 
(see ||82j for an explicit construction of these states). In the large block limit, non- 
zero eigenvalues of the density matrix pt become the same and the density matrix is 
proportional to a 4-dimensional projector 



lim Pb 



Po 



(39) 



Here P4 is the projector onto the 4-dimensional ground space. Both the von Neumann 
entropy and Renyi entropy are equal to In 4 in the limit. 

For finite block, p;,|G; a) = Aq,|G; a), a = 0, 1, 2, 3. The eigenvalues of the density 



matrix are |T9|, Il9 



A. 



3(- 



a 



0; 



(40) 



a 



1,2,3. 



Here A*";, is the size of the block. The matrix A in (l35ll becomes 

/ Ao \ 

Ai 

Ai 

V Ai / 

So that the density matrix pb is the projector P4 multiplied by this matrix A. 



(41) 



6.2. 1- dimensional Homogeneous Model with Generic Spin 

For the high spin homogeneous model in 1-dimension, we have spin-S* in the bulk and 
spin- 1 at both ends to ensure uniqueness of the ground state. The block Hamiltonian 
is 

(Nt-l) terms / 2S \ 

Hb= E E Ajil,l + l)7rjil,l + l)\. (42) 

(l, l+l) g block \J=S'+1 / 

Here projector nj{l, Z + 1) is defined in the same way as in ([2]). Only nearest neighbor / 
and {I + 1) interact, and Aj{l, Z -|- l)'s are positive coefficients. 

As shown in |H2], the ground space of Hb is {S+ 1) ^-dimensional and can be spanned 
by {|VBSArj^(J, M)), J = 0, 1, . . . , S", M = —J, — J+1, . . . , J}. These are also eigenstates 
of Pb with non-zero eigenvalues (see [HS] for an explicit construction of these degenerate 
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VBS states of the block). In the large block limit, all non-zero eigenvalues approach the 
same value and we have 

^m^p,^p^ = ^^^P(5+i).. (43) 

Here P(5+i)2 is the projector onto the degenerate ground space of the block Hamiltonian. 
Both the von Neumann entropy and Renyi entropy are saturated and equal to ln[(5'+l)^] 
in the limit. 

For finite block, pb|VBS7v,( J, M)) = A(J)|VBS7Vi,(J, M)) with eigenvalues of the 
density matrix A (J) independent of M and given by 

m = 1 1 + E(2^ + l)A^-^(/, S)h [IaJ + 1) - Isi^S + 1)) I . (44) 

in which Ii is an Z*^ order polynomial given by a recurrence relation [28], [19], [82], and 
X{1, S) is given by 

{s-ms + i + iy: ^^^^ 

So that matrix A consists of these eigenvalues (I44p in diagonal form and the density 
matrix pb is the projector P(5_|_i)2 multiphed by this matrix A. 

6.3. 1- dimensional Generalized (Inhomogeneous) Model 

For the generalized model in 1-dimension, we label the left ending site of the block by 
/ = 1 with spin value Si = |(Moi + M12) and the right ending site by / = L with spin 
value Sl = \{Ml-i,l + Ml.^ L+i). The block Hamiltonian is 

{Ni-l) terms / Si+Si+i \ 

H,= E Aj(/,/ + l)7rj(/,/ + l) (46) 

(I, l+l) G block \J=Si+Si+-^-Mi^i+i+l j 

with 7rj(/, / + 1) defined in ([2]) and Aj(l, I + 1) positive coefficients. 

As shown in [83], the ground space of Hb is (Mqi + 1)(Ml^l+i + l)-dimensional and 
canbespannedby {IVBSatJJ, M)), J = | J_|, | J_| + 1, . . . , J+, M = -J,-J+1,...,J} 
(see [S3] for an explicit construction of these degenerate VBS states of the block). Here 
J- = |(Moi — Ml,l+i) and J+ = |(Moi + Ml^l+i)- These states are eigenstates of pb 
with non-zero eigenvalues. 

In the large block limit, assuming that Mqi S^ and Mj^ ^.+1 ^ 5"+, we have 

1 

A^r-^oo^"-'^"" {S- + 1YS+ + 1) 

Here S- and S+ are the first and last spins in the block, respectively. The von Neumann 
entropy is equal to the Renyi entropy 

S,.N = Sr = ln[(5_ + 1){S+ + 1)] (48) 

in the limit. 



^lim pb = poo = — — r~rT77^ — r~^P{S-+i){s++i)- (47) 
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For finite block, ph\VBSN^{J, M)) = A(J)|VBSAr,(J, M)). The eigenvalues A(J) are 
independent of M and given by 

n A(/,M,,,+o 



A(J) ^ 




Po/2/.(/,J) , (49) 



(Moi + l)(Mi,L+i + l) 

in which M< = Min{Mj j^i, j = 1, . . . , L — 1} being the minimum of the multiplicity 
numbers, Poly. {I, J) is a polynomial of J depending on /, and A(/, Mjj+i) is given by 

X(IM)- + 



So that matrix A consists of these eigenvalues fl49|) in diagonal form. The density matrix 
is the projector P(5_+i)(s'++i) multiplied by the matrix A. 

6.4- Some Open Problems 

One open problem is the calculation of non-zero eigenvalues of the density matrix pb 
for general and more complicated graphs. One should start with the Cayley tree (also 
known as the Bethe tree). We expect that an exact explicit expression for the non-zero 
eigenvalues is possible, because it has no loops. It is also important to calculate non-zero 
eigenvalues of pb for graphs with loops. We expect that in the large block limit, each 
non-zero eigenvalue should approach the same value = and the entanglement 
entropies should be saturated, i.e. Sy,N = Sr = hi D = In (deg.). 



7. Conclusion 



We have studied the entanglement of the AKLT model. We formulate the AKLT model 
on an arbitrary connected graph. The Hamiltonian ([T]), (I9l) is a sum of projectors which 
describe interactions between nearest neighbors. The condition of uniqueness of the 
ground state relates the spin value at each vertex with multiplicity numbers associated 
with edges incident (connected) to the vertex, see (ITT]) . (IT^ . The unique ground 
state is known as the Valence-Bond-Solid state ([8]), f|T7|) . The VBS state can be used 
instead of the cluster state in measurement-based quantum computation, see [TU]. 

To study the entanglement, the graph is divided into two parts: the block and the 
environment. We investigate the density matrix pb of the block and show that it has 
many zero eigenvalues. We describe the subspace (called the support) of eigenvectors 
of Pb with non-zero eigenvalues. We have proved (see Theorem in Section [5]) that 
this subspace is the degenerate ground space of some Hamiltonian, we call it the block 
Hamiltonian (!23l) . 

The entanglement can be measured by the von Neumann entropy or the Renyi 
entropy of the density matrix pb. Most eigenvalues of pb vanish and have no contribution 
to the entanglement entropies. The density matrix takes the form of a projector on the 
ground space multiplied by another matrix (conjectured in [8^ for an arbitrary graph). 

Non-zero eigenvalues of pb have been calculated for a variety of 1-dimensional AKLT 
models [19], [19], [82], [83]. They are given as illustrations. We find that in these cases 
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the support coincide with the ground space, so their dimensions are equal D = deg. In 
the large block limit, all non-zero eigenvalues become the same and the density matrix 
is proportional to a projector ( l39l) . ( H3i) . ( 1471) . The von Neumann entropy equals the 
Renyi entropy and both take the saturated value 3^,^ = Sr = In D = \n{deg.). 

For more complicated graphs, non-zero eigenvalues of the density matrix are still 
unknown. One open problem is to calculate those eigenvalues. One may start with the 
Cayley tree because there is no loops and we expect to obtain exact exphcit expressions 
of the eigenvalues. 
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Appendix A. Coherent State Basis 

Given the VBS state in the Schwinger representation (|T7j) . it is possible to use a coherent 
state basis [5], [53], [28], [19], [HS]- Then spin operators behave like classical unit 
vectors. The coherent state basis may simplify the calculation of correlation functions, 
eigenvalues and eigenstates of the density matrix. 
We introduce spinor coordinates 

Q \ 

cos-e'2, sin-e-^2 j , < < n, < < 2^. (A.l) 

Then for a point Cl = (sin 9 cos (p, sin 9 sin </>, cos 9) on the unit sphere, the spin-iS coherent 
state is defined as 

\n) = ^ ^ I vac). (A.2) 

Here we have fixed the overall phase (a U{1) gauge degree of freedom) since it has no 
physical content. Note that (lA.2p is covariant under SU{2) transforms (see [22], [52]). 
The set of coherent states is complete (but not orthogonal) such that [1], [28] 

2S+1 ^ ^ 



An 



JdQ\Q){Q\= \S,m){S,m\=hs+i, (A.3) 

m=—S 

where \S, m) denote the eigenstate of S"^ and Sz, and /25+1 is the identity of the (25+1)- 
dimensional Hilbert space for spin-S". The completeness relation (1A.3P can be used in 
taking trace of an arbitrary operator. 

Take the partition function $ = Tr[ |VBS)(VBS| ] for example. Using the coherent 
state basis ( 1A.2I) and realizing that 



(vac|a*+™6^-'^|fi) = J {2S)\ (A.4) 
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we have 



$ = Tr[ |VBS)(VBS| 
(25/ + 1)!" 



n 

. I 



An 



n 

{kl) 



1 — fifc ■ fi; 



(A.5) 



Here (kl) runs over all edges and we are integrating the classical solid angle Cli over each 
vertex /. 



Appendix B. Ground State Degeneracy of the Block Hamiltonian 



We prove Katsura's formula ( l24l) for the degeneracy of ground states of the block 
Hamiltonian. Our proof is based on [48]. The block Hamiltonian is defined in fl23|) . 
We first look at the uniqueness condition (ITTil . For an arbitrary vertex / in the block B, 
the condition can be written as 

2Si = J2Mki = J2MM+J2 Mki, I G B. (B.l) 

k kGB kedE 

Note that the sum over vertices k G dE is outside the block B. These terms are only 
present for boundary vertices / G dB. Expression ( IB. II) is valid for any vertex in the 
block (for a bulk vertex the last summation vanishes). Next we define the block VBS 
state 

\VBSn,)= n (4^1 -^W)^' Vac), keB, leB. (B.2) 

{kl)(^B 

Here edge (kl) lies completely inside the block B. Now an arbitrary ground state of the 
block Hamiltonian Ht takes the following form 

L terms 

n fialbl) IVBS^J, (B.3) 

le dB 

where f{aj, bj) is a polynomial (may depend on the vertex /) in aj and b] and the product 
runs over all boundary vertices (with the number denoted by L). The degree of this 
polynomial is equal to J2kedE ^ki- (Each term in the polynomial must have the same 
total power J2k(^dE ^ki of a] and bj.) It is straightforward to verify that |G) in (IB. 30 is 
a ground state: 

(i) The power of aj and b] in IVBSatj^) is Y^^eB ^ki (see (IB.2P ) so that the total power 



|G) 



of a] and b] in (IB.3P is Y^k&B ^ki + J2k&dE ^ki = 25"; according to (IB.ip . Therefore 



we have the correct power 2*5; of the bosonic operators a\ and b\ for each vertex / 
in the block B (constraint (fT6!) is satisfied); 

(ii) There is no projection on any edge spin value greater than Sk + Si — Mm + 1 because 
of the construction of the block VBS state (\B.2\i . (One could use the same reasoning 
as in Section [3]). 

Therefore the degeneracy deg. of the ground states of Hf, is equal to the number of 
linearly independent states of the form (lB.3p . Since a|'s and 6|'s are bosonic and 
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commute, the number of linearly independent polynomials /(a|,aj) for an arbitrary 
/ is equal to its degree plus one, i.e. {J^k&dE^ki) + 1, V / G dB. So that the total 
number of linearly independent polynomials of the form n/e^a™^ /('^L is the product 
of these numbers for each / G dB. Finally, the ground state degeneracy of the block 
Hamiltonian Hb is (Katsura's formula [|48]) 



deg. 



n 

l&dB 



ykedE 



+ 1 



(B.4) 



In the case of the basic model all M^i = 1, formula flB.4p has a graphical illustration, 
see Figure ([2]). We count the number # of all cut edges (dashed lines) incident to one 
boundary vertex of the block, then add one to the number The degeneracy is the 
product of these (# + l)'s for each boundary vertex. 
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